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The Fermi constant, Gf, describes the strength of the weak force and is determined most 
precisely from the mean life of the positive muon, t^. Advances in theory have reduced the 
theoretical uncertainty on Gf as calculated from to a few tenths of a part per million 
(ppm). The remaining uncertainty on Gf is entirely experimental, and is dominated by the 
uncertainty on Tf^. The MuLan experiment is designed to measure the muon lifetime to part- 
per-million precision, a better-than twenty-fold improvement over the previous generation of 
experiments. In 2007, we reported an intermediate result, = 2.197013(24) fis (11 ppm), 
which is in excellent agreement with the previous world average. This mean life was mea- 
sured using a pulsed surface muon beam stopped in a ferromagnetic target, surrounded by 
a symmetric scintillator detector array. Since this intermediate measurement, the detector 
was instrumented with waveform digitizers, the muon beam rate and beam extinction were 
increased, and two data sets were acquired on different targets, each containing over 10^^ 
muon decays. These data will lead to a new determination of Gf to better than a part per 
million. 



The Standard Model (SM) has provided an excellent description of fundamental particles 
and their interactions since its introduction in the early 1970's. The SM description of the 
weak force requires three input parameters; the most precisely determined are the fine structure 
constant, a, the Fermi constant, Gp, and the mass of the Z-boson, Mz- The Fermi constant, 
Gp, is determined most precisely from a measurement of the positive muon lifetime, r^u, via the 
formula 

where Ag includes corrections from phase space, QED, hadronic and tau loops. In this param- 
eterization, weak and non-weak effects factorize and all weak corrections are included in Gi? as 
Ar, 

where gw is the electroweak gauge coupling. The calculation of the second-order QED cor- 
rections to the determination of Gi^ from reduce the theoretical uncertainty on Aq from 
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~ 30 ppm to less than 0.3 ppm making the uncertainty on the dominant uncertainty on 
Gp- This theoretical progress motivated the current generation of muon lifetime experiments. 

The Muon Lifetime Analysis (MuLan) experiment is designed to measure the muon lifetime 
to part-per- million (ppm) precision using a time-structured muon beam and a symmetric, seg- 
mented detector surrounding a central muon-stopping target. We presented an intermediate 
result, Tn = 2.197013(24) //s (11 ppm), based on data collected in 2004 [2]. Since that measure- 
ment, the experimental electronics were upgraded from discriminators to waveform digitizers, 
and two high-statistics datasets were collected on different stopping targets in 2006 and 2007, 
respectively. 

The experiment runs in a multi-muon mode where several muon decays are measured simul- 
taneously. Muons are provided by a 10 MHz dc beam modulated by a fast-switching, 25-kV elec- 
trostatic kicker [3]. The time-structure imposed by the kicker has a 5-fis "accumulation period," 
during which muons are accumulated in a stopping-target, followed by a "measurement 
period" to observe the decay of the accumulated muons. During the kicker-on measurement 
period, the beam rate is reduced by a factor of ~ 1000. Approximately 50 muons arrive and stop 
in the target during the accumulation period, and 20 survive until the beginning of the mea- 
surement period. An instability in the kicker voltage could affect the observed muon lifetime, 
but the uncertainty on from kicker instabilities has been determined to be less than 0.2 ppm 
for the 2006 dataset and less than 0.07 ppm for the 2007 dataset. 

The positive muons accepted by the beam optics are produced from pion decay at rest and 
have negative helicity. This helicity is preserved as the muons are transported along the beamline 
and stop in the target. A slow precession or relaxation of the average muon spin, combined with 
a possible asymmetry in detector acceptance, could introduce a perturbation on the observed 
muon lifetime. To prevent an unobserved slow spin precession or relaxation from perturbing the 
observed muon lifetime, two targets were chosen to dephase or visibly precess the average muon 
spin. In 2006 a ferromagnetic target with ~ 0.5 T internal magnetic field, called Arnokrome-3 
(AK3) p], dephased the average muon spin during the accumulation period. In 2007, a disk 
of quartz crystal caused ~ 90% of the muons to form muonium, a hydrogen-like bound state 
composed of one muon and one electron. In the spin-0 state, the spin of the muon and the spin of 
the electron are rapidly exchanged, effectively nulling the muon spin effects. In the spin-1 state, 
the spins of the electron and muon are locked together by the hyperfine interaction, and precess 
~ 103 X faster than a bare muon [sj. An array of permanent magnets provide a 130-Gauss field 
to visibly precess the remaining muons. The magnetic field in both target configurations was 
oriented transverse to the beam and initial muon-spin direction. 

Muon decay /i"*" — )• e"*" occurs with ~ 100% probability. Positrons from muon decay in 

the stopping target are observed by the MuLan detector. The detector is composed of a truncated 
icosahedron (soccer ball) configuration of triangular scintillator tile pairs which symmetrically 
surrounds the target (Fig. [T]). This detector segmentation results in low individual tile pair 
count rate, ~ 0.1 hit per tile pair per beam cycle. The detector configuration contains 20 
hexagonal detector housings, each containing 6 tile pairs, and 10 pentagonal detector housings, 
each containing 5 tile pairs. Two pentagonal faces are open, allowing the vacuum pipe to pass 
through the detector. The stopping target is suspended in the center of the detector, inside 
the vacuum pipe. The symmetry of the detector largely cancels any asymmetries from residual 
average muon spin precession or relaxation. 

Several systems are required to read out the scintillation light produced in the tile pairs by a 
through-going positron from muon decay. Each tile is connected to an adiabatic lightguide and 
phototube, and the analog pulses from the phototubes are digitized by 8-bit waveform digitizers 
(WFDs). The WFDs trigger when the analog input rises above threshold, and write out 24 
samples (~ 53 ns) of waveform. Each WFD board has four analog inputs, assigned to a tile 
pair and its symmetrically opposite pair. The 85 WFD boards required for detector readout are 
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Figure 1: Illustrations of the detector, (a) A through-going positron will pass through the inner 
and outer triangular scintillator in a tile pair, (b) A cartoon shows the location of the target 
inside the vacuum beampipe and centered in the detector. 



distributed in 6 VME crates, and each crate is readout by a dedicated frontend computer. 

A backend computer collects the data and controls the time-structure of the data acquisition 
(DAQ). The kicker runs continuously, but the DAQ acquires 5000 beam-cycles at a time, called 
a segment. After each segment, the computers read out how much data is available in the WFD 
FIFO memory, and then triggers the next segment. Data from the previous segment is readout, 
compressed, and assembled by the backend while the current segment is being acquired. This 
parallel structure of the DAQ allows for data rates up to ~ 40 MB/s with no decrease in DAQ 
live-time. The data acquisition is discussed in more detail elsewhere j6|. 

The digitization frequency for the Waveform Digitizers is provided by a master clock, set to 
a frequency of ~ 451 MHz. During data collection and analysis, the precise clock frequency was 
concealed to allow for a blind analysis. During the analysis, times were reported in units of clock 
ticks (ct), where 1 ct~2.2 ns. After the analysis was complete, the true clock frequency was 
revealed. Different blinding offsets were used in 2006 and 2007, and the datasets were analyzed 
independently. 

Overall, ~ 10^^ muon decays were collected on each of the ferromagnetic AK3 and quartz 
targets, respectively. The analysis is performed in two stages. In the first stage, each waveform 
is fit using pulse templates to determine the pulse time and height, as shown in figure [2} In 
the second stage, coincidences are formed between inner and outer tiles of each pair, and the 
coincidences are histogrammed vs. time in measurement period. The histograms are corrected 
for pileup and for multiple hits from single-source events. The histograms are then fit with the 
function 

f{t) = {l + AS{t))Noe-'/-^ + B, (3) 

where AS{t) is small correction for a WFD threshold oscillation, originating in the VME logic 
fanout, A^o is the number of muon decays at t = 0, and B is the fiat background. The magnitude 
of the electronics oscillation correction on the lifetime is 0.60 ppm, and its uncertainty 0.26 ppm. 
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Figure 2: Examples of pulses show the resolving power of the fitting code. When triggered, the 
WFD records 24 samples (~ 53 ns) of waveform data. The raw data is shown in red and the 
fitted pulse template is shown in blue, (a) A fit to a normal single-pulse waveform. Over 90% 
of the data is comprised of similar single-pulse waveforms, (b) The fitter is able to resolve two 
pulses if they are separated by 3 or more samples. 



The lifetime histograms for the two datasets, with fit residuals, are shown in figure [3} 
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Figure 3: Lifetime histograms and fit residuals for the 2006 and 2007 datasets. 



Although the WFDs digitize continually, particles passing through a detector close together 
in time will be incorrectly reconstructed by the analysis, resulting in missed events. This loss of 
hits is called pileup. The pileup is statistically corrected using a shadow-window technique. To 
test the fidelity of the corrections, a set of software deadtimes (ADTs) in the range 5-68 samples 
(11-151 ns) are imposed after a hit at time in beam cycle j. The pileup is corrected by 
searching the same deadtime interval ti to ti + ADT in beam cycle j + 1. If a "shadow" hit is 
found in this interval, it is added back into the lifetime histogram. Several different classes of 
pileup were corrected. A Monte-Carlo simulation of detector hits shows no dependence between 

and deadtime, but a dependency of 0.008 ppm/ns deadtime is observed in the data (Fig. 
|4]). Several effects that could pull vs. deadtime have been ruled-out, and the dependency 
likely corresponds to ~ 0.1% under-correction in the leading-order pileup term. Under-corrected 
leading-order pileup introduces a linear dependence between and ADT. In this case, a linear 
extrapolation to zero deadtime gives the correct value for r^. Overall, the uncertainty on the 
pileup corrections and the extrapolation is 0.2 ppm. 
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Figure 4: Shift in pileup-reconstructed lifetime vs. software-imposed deadtime for the 2006 
dataset. The uncertainty on each point represents the uncertainty of the statistical pileup 
correction. The overall statistical uncertainty of the 2006 dataset is 1.18 ppm. A small slope of 
0.008 ppm per ns deadtime is observed. 

Table 1: Systematic Uncertainties, in units of parts-per- million. These uncertainties represent 
preliminary upper bounds on the uncertainties, and will decrease as studies are finalized. Some 
systematic uncertainties vary between datasets. Correlated uncertainties are denoted by a single 
number, and uncorrelated uncertainties are given with one number for each running year. Effects 
vs. time are evaluated vs. time in measurement period, and effects vs. At are evaluated vs. 
time after a prior pulse. 



Effect 


Size (ppm) 




2006 2007 


Kicker stability 


0.22 0.07 


Spin precession 


n/a 0.20 


Clock calibration 


0.03 


Errant muon stops 


0.10 


Gain stability vs. time 


0.70 


Gain stability vs. At 


0.27 


Timing stability vs. time 


0.09 


Timing stability vs. At 


0.08 


Electronics stability vs. time 


0.26 


Pileup correction 


0.20 


Total systematic 


0.85 


Statistical uncertainty 


1.18 1.7 


Total uncertainty 


1.3 



Other systematic effects are shown in table [T] The largest uncertainty at 0.7 ppm is the 
effect on from the detector response, or gain, vs. measurement time. This effect is still under 



investigation, and the value reported here is an upper hmit. 

Overall, the MuLan experiment has measured to ppm-level precision, and will soon publish 
the final results. Prior to fully unblinding, a relative unblinding, which mapped the two datasets 
into a common blinded space, showed the lifetime difference between the two datasets is 0.3 ppm, 
an excellent agreement. Currently, a few systematic uncertainties are undergoing final analysis, 
and the publication is in preparation. A plot showing the historical precision of Gp and is 
shown in figure [5| including two points in 2010 representing the results from our two datasets. 
After this measurement, the value for will have a precision of 1.3 ppm or better, and Gp will 
have a precision of 0.8 ppm or better. 




Figure 5: Theoretical and experimental uncertainty on Gp vs. time. The second-order quantum 
electrodynamics (QED) corrections to Gp reduced the theoretical uncertainty on Gp from 30 
ppm to 0.3 ppm, making experimental uncertainty on the muon lifetime the dominant contri- 
bution to the uncertainty on Gp. 
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